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Abstract 

The interior of cells is crowded thus making it important to assess the effects of macro- 
molecules on the folding of proteins. Using the Self-Organized Polymer (SOP) model, which 
is a coarse-grained representation of polypeptide chains, we probe the mechanical stability of 
Ubiquitin (Ub) monomers and trimers ((Ub)3) in the presence of monodisperse spherical crowd- 
ing agents. Crowding increases the volume fraction ($c)-dependent average force ((/u(^c)))) 
relative to the value at #c = 0, needed to unfold Ub and the polyprotein. For a given #c, 
the values of (fui^c)) increase as the diameter (cJc) of the crowding particles decreases. The 
average unfolding force (/u($c)) depends on the ratio where D « crdf^)^ with Rg being 
the radius of gyration of Ub (or (Ub)3) in the unfolded state. Examination of the unfolding 
pathways shows that, relative to $c = Oj crowding promotes reassociation of ruptured sec- 
ondary structural elements. Both the nature of the unfolding pathways and (/u($c)) for (Ub)3 
are altered in the presence of crowding particles with the effect being most dramatic for the 
subunit that unfolds last. We predict, based on SOP simulations and theoretical arguments, 
* Corresponding author phone: 301-405-4803; fax: 301-314-9404; thirum@umd.edu 
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that (/u($c)) ~ ^c" , where v is the Flory exponent that describes the unfolded (random coil) 
state of the protein. 

Introduction. 

Cells exist in a crowded environment consisting of macromolecules (lipids, mRNA, ribosome, sug- 
ars, etc.), making it critical to investigate protein folding in the presence of crowding agents.^ 
If the interactions between the crowding agents and the protein of interest are short-ranged and 
non-specific (as is often the case), then the volume excluded by the crowding agents prevents the 
polypeptide from sampling extended conformations. As a consequence, the entropy of the dena- 
tured state ensemble (DSE) decreases relative to the case when the crowding agents are absent. 
These arguments suggest that excluded volume of crowding agents should enhance the stabihty of 
the folded state provided that the crowding-induced changes in the native state are negligible. ^'^ 
The entropic stabilization mechanism, described above, has been used in several theoretical models 
to quantitatively describe the extent of folded protein as a function of the volume fraction, $c, 
of the crowding agents. ^'^ More recently, a theory whose origins can be traced to the concept of 
intra-protein attraction due to depletion of crowding agents near the protein, ^"^ predicts that the 
enhancement in stability, AT($c) = Tf{^c) — Tf{^c = 0) --^ where T/-($c) is the folding tem- 
perature at $c and a is related to the Flory exponent that characterizes the size of the protein in 
the DSE.^ From this prediction it follows that crowding affects the DSE to a greater extent than 
the folded state. Although the precise theoretical predictions of the power law change in AT($c) ^ 
$c changes have not been verified, several experiments using a number of proteins have confirmed 
that indeed T/($c) increases with $c-^^° It cannot be emphasized enough that the theory described 
here applies only to cases when the crowding interactions between crowding agents and proteins and 
between crowding particles themselves are purely repulsive. 

While much less is known about the effects of crowding on the folding kinetics, Cheung et 
al.^ predicted that the entropic stabilization also suggests that the folding rates should increase at 
moderate values of $c- They suggest that crowding can enhance folding rates by a factor e^^^^"^/^'^ 
where A5'($c) ("^ ) decrease in the entropy of the DSE relative to its value in the bulk. From 
the arguments of Cheung et al.^ it follows that the equihbrium changes in the entropy (A5'($c)) of 
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the DSE, with respect to the bulk, should also determine rate enhancement provided that neither 
the barriers to folding^^ nor the native state is perturbed significantly by crowding particles. 

Single molecule force spectroscopy, such as Atomic Force Microscopy (AFM) and Laser Optical 
Tweezers, have been used to monitor the behavior of biopolymers under tension are ideally suited 
to probe the enhancement in crowding-induced stability by a direct measurement of /u($c)- Indeed, 
Ping et al.^^ have recently investigated the effect of Dextran molecules on the mechanical stability 
of (Ub)8. The 8 Ub (Fig. [l|A) modules were N-C linked (i.e., modules i and i+1 were chemically 
linked together in a head-to-tail manner). They found that the average force required to unfold 
a module, (/„($c)), increased by 21% as the Dextran concentration, p, was increased from to 
300 g/L at Vf = 4.2 x 10^ pN/s. Similar results have been obtained recently by Yuan et al. at 
rf = 12.5 X 10^ pN/s.i^ 

Motivated in part by experiments, ^^'^^ we used simulations to investigate the effects of crowding 
agents on the mechanical stability of a protein subject to external tension. We focused on ubiquitin 
(Ub), a 76- residue protein composed of 5 /3-strands and 2 a-helices (Fig. [l]A.), and confined our 
investigation to non-equilibrium 'force-ramp' experiments.^'^ The primary data recorded during 
such an experiment is a trace of the force exerted on the tip as a function of the extension of 
the molecule; a force-extension curve (FEC). When the force exceeds some critical value, the FEC 
displays a sudden increase in length and is often accompanied by a concomitant sharp decrease in 
force. Presumably, the sharp change corresponds to the unfolding of the protein. Typical AFM 
experiments use tandem arrays of proteins which are chemically linked together (often through 
genetic engineering). We use the term module to denote a protein of the array. The FEC resulting 
from such an experiment reveals several equally spaced peaks punctuated by sharp increases in 
the extension of the molecule corresponding to the unfolding of individual modules. The height 
of these force peaks and their shape depend on the loading rate, Vf = kg x v, where kg is the 
cantilever's spring constant and v is the (constant) speed at which the stage is retracted away from 
the cantilever. 

In order to compare to experiments our simulations are performed using coarse-grained models 
for which simulations can be done at rj that are comparable to those in AFM experiments. Our 
work has led to a number of testable results: (1) At $c = 0.3 the average unfolding force for Ub 
increases by at most only 7% compared to $c = 0. We find that (/u($c)) in small crowding agents 
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is greater than in larger particles. (2) In the presence of crowding agents, secondary structural 
elements reform multiple times even after initial rupture. (3) Although large crowding particles 
are predicted to have a smaller effect on (/m($c)) (for a given $c), they can profoundly affect the 
unfolding of poly Ub. We predict that (/u($c)) for a given subunit depends on the number of 
already unfolded portions of the poly protein. This result is important because many naturally 
occurring proteins that are subject to tensile stresses exist as tandem arrays of modules. It further 
suggests that the existence of such redundancy can more properly be understood in the context of 
a crowded cellular milieu. 



Methods. 

Self- Organized Polymer Model for Ub. 

We used a coarse-grained model for proteins to investigate crowding effects on the mechanical 
stability of Ub and (Ub)3 at loading rates that are comparable to those used in AFM experiments. ^^'^^ 
We assumed that Ub could be described using the Self- Organized Polymer (SOP) model; a model 
that has been successfully used to make a number of predictions regarding the unfolding of proteins 
and RNA/®'^^ allosteric transitions in enzymes/^'^^ and movement of molecular motors on polar 
tracks.^'' Previous studies^^ have used more standard Go-models^^'^^ to probe various aspects of 
forced unfolding of Ub. The SOP energy function (Ep) for a protein with N amino-acids, specified 
in terms of the Cq coordinates {i = 1,2, N), is 



Ep = EpENE + -E^" + -^r^ 
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kcal/(mol-nm^), Eh = 1.4 kcal/mol, Ei = 1.0 kcal/mol, and a = 0.38 nm. Note that fc^T ^ 0.6 
kcal/mol ^ 4.2 pN-nm. In Eq. Q Aij = 1 if r°. < 0.8 nm, and = otherwise. Native coordinates 
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r^l is the value of r,-,- in the native structure, k = 2 x 10 
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corresponded to those of the atoms of the 1.8 A resolution Protein Data Bank crystal structure 
1UBQ.24 For Ub iV = 76 and = 228 for (Ub)3. The first term in Eq. Q is the FENE potential^^ 
that accounted for chain connectivity. The second (Lennard- Jones) term accounted for the non- 
bonded interactions that stabilize the native state, and the final (soft-sphere) term accounted for 
excluded- volume interactions (including those of an angular nature). The SOP model is different 
from the Go-model because there are no angular terms in SOP, and the connectivity is enforced 
differently as well. The SOP representation of the polypeptide chain is in the same spirit as other 
coarse-grained models used in polymers. 

Crowding Particles and Interactions with Ub. 

We assumed that the crowding particles are spherical with diameter ctc. ((Tc = 6.4 nm in some 
simulations, while cTc = 1.0 nm in others.) Crowders interacted amongst themselves and with the 
protein, respectively, via the following LJ potentials: 

£„ = 4.,((^)"-(^)Vi)e(C„-r) (2) 

i5, = 4.,((!^)"-(:^)Vi)e(C.,-r). (3) 

where a^p = (cTa + cr/3)/2, ap = a = 0.38 nm, r^^^ = and r^j^ = 2^/^crcp- The Heaviside 

functions truncate the potentials at their minima and thereby ensured that only the repulsive 
portions of Eqs. ^ and ([s]) contributed to interactions involving crowding agent. 

Mimics of Crowding Using Asakura-Oosawa Theory: Even using a coarse-grained SOP repre- 
sentation of proteins, it is difficult to carry out converged simulations in the presence of crowding 
agents. The reason is that the number of crowding agents can be large. Moreover, the separation 
in the spatial and temporal scales of the protein and the crowding particles has to be carefully 
considered to obtain rehable results. In light of these difficulties, it is of interest to consider the 
effective attraction between the sites on the protein using the implicit pairwise potential computed 
by Asakura and Oosawa. The intramolecular attraction arises due to the depletion of crowding 
particles near the protein. To probe the efficacy of these models we employed in some simulations 
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the Asakura-Oosawa model^ of crowding effects. For these simulations, we added the following 
term to the bare SOP Hamiltonian (Eqn. ([T])): 

^AO (r.,) = -^.ksT ^ (^^±^\ ' (i _ + ) a < r < a + a., (4) 

where a = 0.38 nm, = 6.4 nm, $c = 0.3, fc^^ — 4.2 pN-nm, and rjj is the distance separating 
protein beads i and j. 

(Ub)3 Intermodule Interactions. 

For simulations involving (Ub)3, residues in different modules interacted via: 

^PP = ^^(^)'> (5) 

where r is the distance separating the two beads. Note that this potential is short-ranged and purely 
repulsive and that it is the same potential used for non-native intra-protein interactions.^^ 

Simulation Details. 

$c = 0: Hundreds of simulations of 5 x 10^ steps 30/is) at T = 300 K were used to generate initial 
structures for use in the pulling simulations. The protein was completely free in solution (i.e., no 
forces were applied to either terminus), and no crowders were present during the equilibrations. The 
N-terminus of the protein was subsequently translated to the origin, and the protein was rotated such 
that its end-to-end vector, R, (i.e., the vector pointing from the N-terminal bead to the C-terminal 
bead) coincided with the pulling (+z) direction. 

An unfolding trajectory was initiated by selecting a random initial structure from amongst the 
set of thermally equilibrated structures, and tethering a harmonic spring to the C-terminal bead. 
The N-terminal bead was held fixed throughout the simulations. Tension was applied to the protein 
by displacing the spring along the +z axis and resulted in application of the following force to the 
C-terminal bead: 

/. = -K {[z{t) - ^(0)] - [z,{t) - zM]) , (6) 
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where kg is the spring constant, z{t) — R(i) -z, and Zs{t) corresponds to the displacement of the end of 
the spring. Note that ks was also used to constrain the simulation to the z-axis; f.^ = —ks[x{t)—x{0)] 
and fy = —ks[y{t) — y{0)]. The displacement of the spring was updated at every timcstcp. 

We simulated forced-unfolding of monomeric Ub at four different r/ (160 x 10^ pN/s, 80 x 10^ 
pN/s, 20 X 10^ pN/s, and 4 x 10^ pN/s), while simulations on N-C-linked (Ub)3 were performed at 
rf — 640 X 10^ pN/s. All overdamped force-ramp simulations were performed at the same speed 
V = 10312 nm/s, and spring constants were varied over a range from 0.3879 pN/nm - 31.032 pN/nm 
to achieve the aforementioned Vf (via the relation r/ = ksv). Our simulations were realistic because 
they maintained loading rates consistent with experiment and because r/ is the prime determinant 
of unfolding pathway. 

$c 7^ 0.0: Simulations involving explicit crowders were carried out at a fixed volume fraction 
$c = 0.3 and with a fixed number Nc = 100 of crowding spheres. {Nc was fixed to render the 
problem computationally tractable). Using the relation $c = adjusted the length of 

a side of the cubic simulation box ( L ) to maintain $c — 0.3. Thus, L — 35.8 nm when ctc = 6.4 
nm and L — 5.6 nm when = 1.0 nm. Explicit crowders were added to the simulation after 
loading an equilibrated structure but before the application of tension. Initial crowder positions 
were chosen randomly and in a serial manner from a uniform distribution. If the distance between 
an initial crowder position and that of another crowder or protein bead did not exceed the sum of 
their radii, then the prospective position was rejected and another random position chosen to avoid 
highly unfavorable steric overlaps. 

Periodic boundary conditions (PBC) and the minimum image convention^^ were employed in 
the simulations. Two sets of coordinates were stored for protein beads at every timestcp; PBC 
were applied to one set and the other was propagated without PBC. Distances between protein 
beads were calculated from the uncorrected set of coordinates without minimum imaging, while 
protein-crowder distances were calculated from the PBC coordinates with minimum imaging. 

To improve simulation efficiency, a cell list^^ was used to calculate crowder-crowder and protein- 
crowder interactions. The entire simulation volume was partitioned into 64 subvolumes, and it was 
only necessary to calculate interactions within a subvolume and between beads of the subvolume and 
those of 13 of its 26 neighbors. The cell-list was updated at every timestep to ensure the accuracy 
of the simulations. 
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The equations of motion in our overdamped simulations (used in all force-ramp simulations) were 
integrated with a timestep h = O.Olri^ {tl = 2.78 ps) using the method of Ermak and McCammon.^^ 
The friction coefficient of the crowders, was determined via the relation ^ = where cxc is the 
crowder diameter, a = 0.38 nm is the diameter of a protein bead, and ( = 83.3 x [^) = 9 x 10~^ g/s 
is the friction coefficient associated with a protein bead of mass m = 3 x 10~^^ g. Simulated-times 
were translated into real-times using th = (^y) '^l'x (^)-^^ At T = 300 K, th = 543.06 ps, and 
since h = 0.01 x tl the real-time per step is 5.4306 ps. 



Results and Discussion. 

Monomeric Ub at $c = 0.0. 

At $c = 0.0, forced unfolding of Ub was simulated at four different r/ (4 x 10^ pN/s, 20 x 10^ pN/s, 
80 X 10^ pN/s, 160 X 10'^ pN/s), where the lowest value corresponds approximately to the value used 
in the pulling experiments of Ping et al.,^^ and all Vf are experimentally accessible. 

Force Profiles: Fig. [2]^ A and B) provides examples of FEC's collected at the highest and lowest 
rf. We used a nominal contour length of (A^ — 1) cr = 75 x 0.38 nm and unfolding forces, /„, to 
determine contour-length increments, AC, for each trajectory at Vf = 160 x 10^ pN/s (Fig. [2|A). We 
identified fu with the peak of the EEC before the stick-slip transition. The average extension 
(A£) = 23.991 ±0.010 nm is in excellent agreement with the experimental result of 24 ± 5 nm found 
by Carrion- Vazquez et al.^'' The projection (z^) of the end-to-end vector at /„ in the z-direction 
varied between 4.1 nm and 4.7 nm, depending on rj. Since the native end-to-end distance, zq = 3.7 
nm, Zu — Zq = Azu ranges from 0.4-1.0 nm. The lower end of this range is slightly larger than 
the 0.25 nm transition-state distance for the mechanical unfolding of the structurally similar titin 
immunoglobulin domains. Indeed, we expect Az^ > 0.25 nm, because of the non-equilibrium 
nature of the simulations. Larger rj typically lead to larger (Azu) (~ ksT/ {fu) ln(r/)). 

Average unfolding forces, {fu{^c)), depended approximately logarithmically on r/^^ (Fig. |2]C), 
and {fui'f'f = 4 X lO'^pN/s)) = 136 pN is in fair agreement with the experimental value of 166 ± 33 
pN observed by Ping et al.^^ at Vf = 4.2 x 10^ pN/s. The (zu) also showed a logarithmic dependence 
on r/, but the difference between the value calculated at r/ = 160 x 10'^ pN/s and that calculated at 
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rj = 4x 10^ pN/s is small (~ 2 A). Although the underlying free-energy landscape is time-dependent 
in a non-equilibrium force-ramp pulling experiment, these results suggest that the distance from the 
native state to the transition state is small. It is likely that at loading rates that are achieved in 
laser optical tweezer experiments (~ 10 pN/s), the location to the transition state would increase 
because the response of biopolymers to loading rate changes from being plastic (low rj) to brittle 
(high Tf).^^ 

Unfolding Pathways: In the dominant pathway, unfolding proceeded in a fairly Markovian fash- 
ion with the primary order of events following the sequence [31/ (35 (33/f35 (33/(34 (31/(32 
(Figs. [T]B and [3]). This is precisely the same sequence seen in the simulations by Li et al.^^ Alter- 
native pathways, reminiscent of kinetic partitioning^^ observed in forced-unfolding of GFP^^ and 
lysozyme,^^ were also infrequently sampled. For example, at rj = 4 x 10'^ pN/s, ~ 6% of the tra- 
jectories unfolded as follows: pi/ (35 (31/(32 (33/(35 (33/(34 (Fig. |i|B), while the remaining 
~ 94% followed the dominant pathway. At the highest loading rate, one frequently observed the 
following sequence of events pl/P5 (33/(35 (31/(32 (31/(32 (reform) ^ (33/(34 (31/(32, 
where (31/(32 ruptured but then reformed prior to the rupture of (33/(34. As illustrated in Fig. [2]^ A 
and B), unfolding events at smaller rj tended to result in larger molecular extensions. 

The non-equilibrium character of a pulling experiment decreased with decreasing r/, and smaller 
Tf resulted in smaller force-drops after the unfolding force, /„, is reached. Since the force applied by 
the spring to the end of the protein did not fall off as sharply at lower rj, more of the protein was 
extended during an unfolding event. The smaller extensions following unfolding events at higher rj 
(relative to those observed at lower rj) were responsible for the (31/(32 unfolding/refolding events 
mentioned above because the applied tension was very low after the initial rupture event (Fig. 
and B)). At lower rf, this situation no longer held because the initial unfolding event resulted in a 
chain extension that was a significant fraction of the chain's contour length (Fig. |2]B). 

Crowding Effects on Ub ($c = 0.3). 

Depletion forces stabilize proteins and shift the folding equilibrium towards more compact states.'^'' 
These forces result from an increase in the entropy of the crowding agents that more than com- 
pensates for an increase in the free-energy of a protein molecule upon compaction. Simulations of 
forced-unfolding of Ub in the presence of explicit crowders of diameters Uc = 6.4 nm and cxc = 1.0 
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nm were used to assess the contribution of the depletion forces to mechanical stability. Sixteen tra- 
jectories were collected for each rj investigated. Three Vf (20 x 10^ pN/s, 80 x 10^ pN/s, 160 x 10'^ 
pN/s), were explored for the o"c = 6.4 nm sized depletants. Only the two highest r/ (80 x 10'^ pN/s, 
160 X 10^ pN/s), were explored for the cXc = 1.0 nm sized crowders. 

Small crowding particles increase unfolding forces: Example FEC's collected in the presence 
of crowders of diameter cTc = 6.4 nm and cTc = 1.0 nm are illustrated in Fig. |4]^A and B). The 
(Tc = 6.4 nm curves in Fig. |4] look qualitatively very similar to those seen at $c = 0.0, while the 
FEC's collected at cTc = 1.0 nm look qualitatively different. For example, larger (fy) than those 
observed at either $c = 0.0 or in the presence of the cTc = 6.4 nm crowders are apparent in the 
FEC's collected at cxc = 1.0 nm (Fig. |4|^A and B)). Indeed, Fig. |4]C reveals that this observation 
is quantitatively accurate. Although the (/„) in the presence of the Uc = 6.4 nm crowders were 
statistically indistinguishable from the (/„) at $c = (compare Fig. |4]C with Fig. [2]C), the average 
unfolding forces in the presence of the Cc = 1.0 nm crowders were statistically greater than those 
measured at $c = 0.0. At Vf = 160 x 10^ pN/s, the (/„($c = 0.3)) in the presence of the cXc = 1.0 
nm crowders exceeded that at $c = 0.0 by 3%, while at = 80 x 10^ pN/s the increase was 4% 
(compare Fig. |4]C with Fig. |2]C). In one respect these results are not surprising; it follows from the 
AO theory and Eq. ([s]) that smaller crowders stabilize Ub more than larger ones. On the other 
hand, the extent of stabilization (as measured by increases in (/«)) was small. We should emphasize 
that although the increase in the unfolding force is small, the stability change upon crowding is 
significant. The enhancement in stability is AG ~ {fu{^c)){zDSE) ~ 5 /e^T using a 3% increase in 
the unfolding force, and {zdse) the location of the unfolded basing ^ 5 nm. 

Crowding leads to transient local refolding: The unfolding pathways in the presence of crowding 
agents of both sizes were very similar to those seen at $c = 0.0. Nevertheless, at $c = 0.3 there 
tended to be more unfolding/refolding events as the molecule extended past Zu (Fig. |3]) than at 
$c = 0.0. As illustrated, strand-pairing between (31 and /52 and between (33 and /34 persisted to 
a greater extent after the initial unfolding event at $c = 0.3 than at $c = 0.0. As Ub passed 
through the point {zu,fu), its termini were extended to distances greater than the diameter of even 
the larger crowders. Depletion forces resulting from the presence of 6.4 nm crowders act on these 
larger length scales. Assuming that the tension applied to the C-terminus is small enough (e.g., 
after rupture events at higher loading rates), then such depletion forces can promote reformation of 
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contacts between secondary structural elements several times during the course of a trajectory. If 
this is indeed the case, then it suggests that unfolding polylJb may be different from the unfolding 
of monomeric Ub, because depletion effects should increase with number of modules in the tandem 
(see below). 

AO Model For Forced- Unfolding in the Presence of Crowders: We used the Asakura-Oosawa AO 
model^"^ (Eq. Q) of the depletion interaction to model the effects of a crowded environment on 
Ub. The AO theory has been successfully used to model the effects of a crowded environment on 
polymers and colloids,^' '''^^'^^ and we found that it gives qualitatively accurate results for Ub. We 
used Eq. ^ to approximate the effective interaction between spherical protein beads immersed in a 
crowded solution of volume fraction $c = 0.3. From the form of the AO potential between protein 
beads, it follows that the range of the potential is proportional to cTc, but the strength is greater 
for smaller crowders. Indeed, as revealed in the previous section, simulations in the presence of 
explicit crowders showed that 1.0 nm crowders resulted in larger average unfolding forces than 6.0 
nm crowders (Fig. |4]C). 

Although the AO-potential yields qualitatively accurate results, use of the AO-potential of Eq. 
^ to implicitly model non-bonded interactions in Ub did not yield quantitatively accurate results. 
At r/ = 160 X 10^ pN/s (/„($c = 0.0)) = 175. 78± 1.52 pN, while {fu{^c = 0.3)) = 266.26±2.70 pN. 
Thus, simulations with the AO potential led to a mean unfolding force that is roughly 50% greater 
than in its absence. This disagrees sharply with the (/u($c = 0.3)) = 173.64 ± 3.49 pN resulting 
from our own simulations in the presence of explicit crowding agent at $c = 0.3 at Vf = 160 x 10^ 
pN/s ( Fig. |4]C ). Indeed, the result also stands in marked contrast to the experimental results 
of Ping et al.^^ on octameric Ub, which saw a maximum increase in (/„) of 21% (at $c > 0.3, 
r/ = 4200 pN/s, and with ~ 7.0 nm) over the (/„) = 166 pN at <l>c = 0.^^ 

There are a couple of origins to the discrepancy. First, the AO-potential was derived to un- 
derstand the equilibrium of colloidal spheres and plates in the presence of smaller-sized spherical 
crowding agents. Our experiments were of a non-equilibrium nature, so it is somewhat unreasonable 
to expect such simulations to yield quantitatively accurate unfolding forces or dynamics. Second, as 
pointed out by Shaw and Thirumalai,^^ three-body terms are required to properly model depletion 
effects even in good solvents let alone in the poor-solvent conditions of our simulations. To elabo- 
rate, let us consider the volume excluded to crowders by a Ub molecule, Vex(Ub), to be the volume 
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enclosed by a union of spheres of radii Si = With an AO potential, the volume excluded to 

the spherical crowding agents is: 



where N is the number of residues in monomeric Ub (76), V{Si) is the volume of Si and V{Si fl Sj) 
is the volume associated with the overlap of Si and 5*-,-. Eq. ([7| neglects the overlap of three 
or more spheres. The importance of such overlaps, for soft-spheres, increases as the crowders 
become much larger than the protein beads, and the neglect of such overlaps is the reason for 
the quantitative inaccuracy. As the size of the crowders increases (i.e., as the thickness of the 
depletion layer surrounding the protein beads increases), the surface enclosing Ub becomes more 
spherical, loses detail, and undoubtedly changes less in response to changes in the conformation 
of the molecule. Since depletion forces are proportional to the change in Vj;(Ub) with respect to 
changes in Ub conformation, Eq. ([t]) overestimates the size of depletion forces in the presence of large 
crowders. Despite these limitations the AO model, which is simple, can be used to provide qualitative 
predictions. Finally, we note that experiments typically use polyproteins to study force induced 
unfolding (e.g.. Ping et al.^^ used (Ub)8 in their experiments). For polyproteins the quantitative 
accuracy of the AO theory for unfolding in the presence of crowders of diameter cTc = 6.4 nm is likely 
to increase, because the volume excluded to the crowders will change significantly with changes in 
the conformation of the polyprotein. 



From arguments based on volume exclusion that lead to crowding-induced entropic stabilization of 
the folded structures, it follows that crowding effects should be more dramatic on poly Ub than the 
monomer. In order to illustrate the effect of crowding on stretching of (Ub)3 we chose = 6.4 nm, 
which had negligible effect on (/„($c)) for the monomer. However, we found significant influence 
of the large crowding particles when (Ub)3 was forced to unfold in their presence. The FEC (Fig. 
[5]) shows three peaks that corresponded to unfolding of the three domains, when simulations were 
performed at crowder volume fractions $c = 0.0 and = 0.3 at rj = 640 x 10^ pN/s. Fig. |6] 
presents average unfolding forces as a function of the unfolding event number. Although the first 
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(Ub)3 at $c = 0.0 and $c = 0.3. 
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two unfolding events were statistically indistinguishable at $c = 0.0 and $c = 0.3, the final event 
occurred at much larger (/„) in the presence of crowders than in their absence. 

Order of Unfolding was Stochastic: (Ub)3 has 3 chemically identical modules. In the pulling 
simulations the N-terminus of module A was held fixed, while force was applied to the C-terminus 
of module C. Figure [7] illustrates the time-dependence of contacts between secondary structure 
elements of modules A, B, C. It is clear from Fig. [7]A. that module C unfolded first, followed 
by module 5, and finally by module A. On the other hand, the order of events in Fig. [7^ was 
A ^ C ^ B. The frequency with which the 3! = 6 possible permutations of these orders at •I'c = 
and at $c = 0.3 were observed is presented in Table [TJ It is clear that at both $c = and at 
$c = 0.3, the most probably order of events was C ^ B —>■ A. At $c = this order was only 
marginally more probable than the order C —>■ A ^ B, while at $c = 0.3 C ^ B —>■ A became 
overwhelmingly more probable than any other unfolding order. In none of the simulations at $c = 
or $c = 0.3, did the rupture of module C occur as the final event. 

Unfolding Within a Module Depended on Proximity to the Point of Force Application: At Vf = 
640 X W^pN/s, (Ub)3 is fairly brittle, and the rupture of contacts within a module occured nearly 
simultaneously. Nevertheless, by carefully examining time-dependent contact maps such as those 
illustrated in Fig. [7| we were able to determine (1) at both $c = and $c = 0.3, /51//55 contacts 
were the first to rupture, and (2) only for module C was this rupture event invariably followed by 
the loss of (33/ (35 contacts. When other modules ruptured, loss of (31/(35 contacts was occasionally 
followed by loss of the (31 / (32 strand-pair contacts. 

(Ub)3 Must Achieve a Larger Rg to Rupture at = 0.3: Figure |8]A. illustrates the time depen- 
dence of {Rg) at $c = and at = 0.3. Interestingly, the plot reveals that after the second rupture 
event, the (i?g($c, t)) increased more rapidly in the presence of crowders than in their absence. This 
is likely a refiection of the fact that at $c = modules A and C were the first two modules to unfold 
in 44% of the trajectories, while at $c = 0.3 these two modules were the first to unfold in only 19% 
of trajectories. Thus, {Rg{^c,t)) increased more rapidly in the presence of crowders, because the 
Rg of (Ub)3 with two adjacent modules unfolded is larger than that with two unfolded but non- 
adjacent modules. Interestingly, these differences are masked in the time-dependent increase of the 
end-to-end distance (Fig. [8^). Figure [sjA. also reveals that the horizontal inflection points marking 
the third unfolding event occur at different times at $c = 0.3 and in the absence of crowders. The 
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difference between these two times was responsible for the difference in average unfolding forces of 
^ 14 pN illustrated in Fig. |6j Thus, it is clear that despite the highly non-equilibrium nature of 
this pulling experiment, depletion effects were substantial and it required a much greater force to 
reach in the presence of these crowders than in their absence. 

We predict that systematic experiments will reveal that polylJb molecules composed of larger 
numbers of modules will show greater increases in {fu{^c)) relative to (/m($c = 0)) than polyUb 
molecules composed of fewer repeats. The size of these differences should increase with decreasing 
loading rate. Finally, it may even be possible to observe differences in the (fu) as a function of 
unfolding event number (as in Fig. [6]). The increase in (/u($c)), at a fixed $c, for poly Ub is likely 
to be even more significant for small crowding agents as shown in Eq. (|8| (see below for further 
discussion) . 

Conclusions. 

General theory, based on the concept of depletion effects (see Eqs. ^ and ([s])), shows that crowding 
should enhance the stability of proteins, and hence should result in higher forces to unfold proteins. 
However, predicting the precise values of {fu{^c)) is difficult because of the interplay of a number of 
factors such as the size of the crowding agents and the number of amino acid residues in the protein. 
Despite the complexity a few qualitative conclusions can be obtained based on the observation that, 
when only excluded volume interactions are relevant, then the protein or polyp rotein would prefer 
to be localized in a region devoid of crowding particles. The size of such a region D ^ cTci-f^)^ ■ 
If D ^ Rg then the crowding would have negligible effect on the unfolding forces. The condition 
D ^ Rg can be realized by using large crowding particles at a fixed $c- In the unfolded state, 
Rg ~ 0.2A^'''^ nm^^ which for Ub leads to Rg ~ 2.7 nm. Thus, = OAa^- These considerations 
suggest that the crowder with cTc = 6.4 nm would have negligible effect on the unfolding force, 
which is in accord with the simulations. On the other hand, -U- ~ 0.4 when ac = 1 nm, and hence 
we expect that the smaller crowders would have measurable effect on the unfolding forces. Our 
simulations are in harmony with this prediction. We expect that for the smaller crowding agent 
ifui^c)) would scale with $c in a manner given by Eq. (|8]). In general, appreciable effect of crowding 
on the unfolding forces can be observed only for large proteins or for polyproteins using relatively 
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small crowding agents. 

Although we have only carried out simulations for Ub and (Ub)3 at one non-zero ^c, theoretical 
arguments can be used to predict the changes in (/„($c)) as $c increases. The expected changes in 
the force required to unfold a protein can be obtained using a generalization of the arguments of 
Cheung et al.^ In the presence of crowding agents the protein is localized in a region that is largely 

1 /3 

devoid of the crowding particles. The most probable size of the region is D ~ o"c$c where cxc 
is the size of the crowding agent. If the structures in the DSE are treated as a polymer with no 
residual structure then the increase in entropy of the DSE upon confinement is AS/ks ~ {Rg/DY^'^ 
where Rg is the dimension of the unfolded state of the protein. If native state stabilization is solely 
due to the entropic stabilization mechanism we expect: 

~ TAS/L, ~ (8) 

where /n($c) is the critical force for unfolding the protein, Lc is the gain in contour length at the 
unfolding transition, and z/(^ 0.588) relates Rg to the number of amino acids through the relation 
Rg ~ o,£,N^ (ajj varies between 2-4 A). A few comments regarding Eq. ^ are in order. (1) The 
$c dependence in Eq. ([s]) does not depend on the nature of the most probable region that is free 
of crowding particles. As long as the confining region, which approximately mimics the excluded 
volume effects of the macromolecule, is characterized by a single length D, we expect Eq. (|8]) to 
be valid. (2) The additional assumption used in Eq. Q is that ^ 1, and hence there may be 
deviations due to finite size effects. (3) The equivalence between crowding and confinement breaks 
down at large $c values. Consequently, we do not expect Eq. (|8| to fit the experimental data at 
all values of $c- (4) It follows from Eq. ([s]) that, for a given Rg, small crowding agents are more 
effective in stabilizing proteins than large ones. Thus, the prediction based on Eq. [8] is supported 
by our simulations. (5) From the variation of {fu{^c)) with $c Ping et al.^^ suggest that (/„) ~ 
However, the large errors in the measurements cannot rule out the theoretical prediction in Eq. ([s]). 
We have successfully fit their experimental results using Eq. ([s]) (Fig. |8]C). Additional quantitative 
experiments are required to validate the theoretical prediction. 

It is difficult to map the concentrations in g/L used in the study of Ping et al.^^ to an effective 
volume fraction because of uncertainties in the molecular weight of Dextran used in the study. 
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Hence, a quantitative comparison between theory and experiments is challenging. A naive estimate 
may be obtained using the values reported by Weiss et al.^^ Given that the Dextran used in the study 
is thought to have an average molecular weight of 40 kDa and an estimated average hydrodynamic 
radius of 3.5 nm,^^ we find that p = 300 g/L corresponds to a volume fraction $c = 0.8, which 
is very large. Nevertheless, $c must be large when p = 300 g/L. Alternatively, we estimated (Tc/2 
for Dextran using crc/2 ~ udN^^^ where is the number of monomers in a 40 kDa Dextran is 
40/0.162 ^ 147. If the monomer size a/? ~ 0.4 — 0.45 nm, then we find $c ~ 0.3 — 0.4. If we assume 
that p = 300 g/L corresponds to ^ 0.4 and that (/n($c)) ~ ^c^^ (Eq. (§)), then at = 0.3 
we would expect a nearly 18% increase in (/u). Similarly, if p = 300 g/L corresponds to $c — 0.3 
then we would expect an increase in (fu) of approximately 21%. In any case, we can say that at 
physiologically relevant volume fractions (<l>c £ [0-1,0.3]), the percent increase in (/„) is likely to be 
< 20%. Our simulations for ac = 1.0 nm predict an increase of 3-4%, which shows that a more 
detailed analysis is required to obtain an accurate value of cTc for Dextran before a quantitative 
comparison with experiments can be made. The larger increase seen in experiments may also be a 
reflection of the use of (Ub)8 rather than a monomer. 

Regardless of the crowder size we find that the unfolding pathways are altered in the presence 
of crowding agents. It is normally assumed that the rupture of secondary structure elements is 
irreversible if the applied force exceeds a threshold value. However, when unfolding experiments are 
carried out in the presence of crowding particles, that effectively localize the protein in a smaller 
region than when $c = 0, reassociation between already ruptured secondary structures is facilitated 
as shown here. Thus, forced-unfolding cannot be described using one dimensional free energy profiles 
with the reaction coordinate. 

We find that the average unfolding force for the final rupture event of the unfolding of (Ub)3 
occurred at much larger values in the presence of crowders than in their absence. With o"c = 6.4 
nm, which has practically no effect on the unfolding force of the monomer, and $c = 0.3 even with 
unfolding of two modules the interactions between the stretched modules and protein are small 
{D ^ 1.2(Tc). Only upon unfolding of the third Ub do crowding effects become relevant, which 
leads to an increase in (/„($£)). Our results suggest that {fui^c)) / {fui^)) should increase with 
the number of modules in the array and that it may be possible to detect differences in (/„) which 
are conditional on the unfolding event number. We speculate that naturally occurring polyproteins 
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that are subject to mechanical stress have evolved to take advantage of precisely such enhanced 
depletion effects. 
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Tables 



Table 1: Module Unfolding Order Frequencies at $c = and $c = 0.3. 



Unfolding Order 


Frequency Observed at $c = 


Frequency Observed at $c = 0.3 


C ^ A 


0.44 


0.56 


C ^ A^ B 


0.38 


0.13 


B^A^C 


0.00 


0.00 


B^C A 


0.13 


0.25 


A^B^C 


0.00 


0.00 


A^C^B 


0.06 


0.06 
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Figure Captions. 

Fig. [l| (A) Cartoon representation of the native structure of ubiquitin (PDB accession id lUBQ) in 
the presence of spherical crowding agents. The five beta-strands, labeled /31 through /35, are colored 
in yellow. The two alpha-helices (al and a2) are shown in purple. The N- and C-terminal beads are 
represented as spheres. In our simulations the N-terminal bead was held fixed while the C-terminal 
bead was pulled via a tethered spring. (B) Snapshots from an unfolding trajectory illustrating the 
main ubiquitin (Ub) unfolding pathway (brown-dashed arrows) and an alternate unfolding pathway 
(green-dotted arrow). In both pathways the initial unfolding event corresponds to separation of 
the C-terminal strand /35 from the N-terminal strand (31. Along the main pathway, this is quickly 
followed by separation of (35 from (33. The penultimate rupture event along the main pathway 
corresponds to disruption of the (33/(34 strand-pair, while the N-terminal (31/(32 strand pair is the 
last to break. The trajectory illustrated here was generated at $c = 0.0 and r/ = 160 x 10^ pN/s. An 
alternate pathway was observed at $c = 0.0 and rj = A x 10^ pN/s. Along this pathway separation 
of (35 from (31 is followed by separation of (31/(32. The final two rupture events correspond to those 
of the (35/(33 contacts and (33/ (3A strand-pair respectively. (Figures generated with VMD^^) 
Fig. [2| Force-extension curves (FEC's) at two different loading rates, (A) rf = 160 x 10^ pN/s, and 
(B) = 4 X 10^ pN/s. Data from the simulation is presented as a red trace. For each trajectory 
a black arrow points to the unfolding force, /„. Zu corresponds to the extension of the molecule 
along the pulling (i.e., z-) direction evaluated at /„. AC (dotted blue line) is the contour length 
increment, and is a measure of the amount of chain released in an unfolding event. We measured 
A£ as £ — Zu, where C = (A — 1) x cr = 75 x 0.38 nm is a nominal contour length of 28.5 nm. Stars in 
each subfigure mark the minimum force observed after an unfolding event, and chain conformations 
corresponding to the starred points in figures A and B are illustrated at the center of the figure. 
Unfolding events at smaller Vf resulted in larger molecular extensions before significant resistance 
was encountered. (Yellow arrows correspond to beta-strands and purple cylinders correspond to 
a-helices). (Figures generated with VMD"^^). (C) (/„) vs. rj evaluated at $c = 0. The red curve 
corresponds to a linear-least squares fit to the set of basis functions {l,ln(rj)} and demonstrates 
that (/„) ~ In(rj). (Note that the abscissa is a log-scale). Each point is labeled Mean ± standard 
error. Statistics at r/ = 160 x 10^ pN/s, 80 x 10^ pN/s, 20 x 10^ pN/s, and 4 x 10^ pN/s were 
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calculated from 50, 49, 50, and 16 trajectories respectively. 

Fig. [sj Rupture events at $c = 0.0 (A) and at $c = 0.3 and Gc = 6.4 nm (B). The figure illustrates 
that after the initial rupture event at step 0, subsequent unfolding was affected by the crowding 
agent. (31/(32 contacts and /53//34 persisted longer and there were many more local refolding events 
in the presence of the crowders (B) than in their absence (A). (C) Snapshots from an unfolding 
trajectory at $c = 0.3 and r/ = 4 x 10^ pN/s illustrate the primary difference between unfolding 
at $c = 0.3 and at $c = 0.0. The brown-dotted arrow shows unfolding without any local-refolding 
events, while the sequence of black arrows illustrates local-rupture and -refolding events. The hall- 
mark of forced-unfolding in a crowded environment is repeated breaking and reforming of contacts 
after an initial rupture event. (Figures generated with VMD^^) 

Fig. |4j (A) Examples of force-extension traces resulting from simulation in the presence of spherical 
crowding agents of diameter o"c = 6.4 nm and obtained at rj = 80 x 10^ pN/s. (B) Examples of force- 
extension traces resulting from simulation in the presence of spherical crowding agents of diameter 
(Tc = 1.0 nm and obtained at = 80 x 10^ pN/s. Both subfigures are labeled as in Fig. [2]^A and 
B). (C) {fu) vs. Tf evaluated at $c = 0.3. Black triangles and red circles correspond to spherical 
crowders of diameter ctc = 1.0 nm and Uc = 6.4 nm respectively. Each point is labeled Mean ± 
standard error. Statistics for each point were calculated from 16 independent trajectories. Only the 
o"c = 1-0 nm had an appreciable effect on (/„) when compared to those obtained at identical and 
at <l>c = 0.0 (see Fig. |2|C) 

Fig. |5| EEC's for (Ub)3 forced unfolding at $c = 0.0 (A) and at $c = 0.3 (B). Trajectories were 
generated at rj = 640 x 10^ pN/s and with crowders of diameter cXc = 6.4 nm. Black arrows mark 
each trajectory's three unfolding events. 

Fig. |6| (/„) vs. unfolding event number for the unfolding of (Ub)3 in the presence of spherical 
crowders of diameter ac = 6.4 nm ($c = 0.3, black triangles) and in their absence ($c = 0.0, 
red circles). The individual modules of the polylJb tandem were N-C-linked and the loading rate 
was 640 X 10"^ pN/s. Each point is labeled Mean ± standard error. Statistics for each point were 
calculated from 16 independent unfolding trajectories. An unfolding event corresponded to the 
unfolding of an individual module. Note that although the crowders had little effect on (/„) for the 
first and second unfolding events, they had a substantial effect on the last unfolding event. (/„) 
increased by ^ 14 pN for the last unfolding event. 
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Fig. [Tj Illustration of the stochastic nature of module unfolding. In (A) module C (the most 
proximal to the applied force) unfolds first, followed by module 5, and finally by rupture of module 
A. In (B) the order of module unfolding events is A ^ C ^ B. The two trajectories illustrated 
here were both collected at $c = and rj = 640 x 10^ pN/s. Table [l] provides the frequencies at 
which the different possible orders were observed at $c = and at $c = 0.3. The diameter of the 
crowders is dc = 6.4 nm. 

Fig. [8] (A) {Rg{t)) versus time at <l>c = (red) and at <l>c = 0.3 (black). Although the (Rg) 
increases more rapidly with time after the second rupture event at = 0.3 than at = 0, the 
inset reveals that a larger Rg must be achieved to initiate the final rupture event in the presence of 
crowding particles with (Tc = 6.4 nm. (See text for additional discussion). (B) {z{t)) versus time at 
$c = (red) and at $c = 0.3 (black). z{t) cannot discriminate between unfolding at the different 
volume fractions and hence is less suitable (than Rg) as a potential reaction coordinate. (C) The 
experimental results of Ping et al.^^ (red circles) for unfolding-force, (/„) as a function of Dextran 
concentration, p. The black line is a fit assuming (/„) ~ p and the blue assuming (/„) ~ p^/^. 
Although both fits are consistent with the data, based on theoretical considerations (see Eq. (|8])) 
we prefer the blue fit (see text). (Standard deviations taken from Ping et al.^^) 
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Figure 3: 
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